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PHOTOS as a pocket parton shower: 
flexibility tests for the algorithm 
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^ I Recently version 2.14 of the PHOTOS Monte Carlo algorithm, written for bremsstrahlung 

(y-^ I generation in decays became available. In ref. Q, detailed instructions on how to use the 

^ ■ program are given. With respect to older versions |l2l|3l of PHOTOS, it now features: improved 

^ ■ implementation of QED interference and multiple-photon radiation. The numerical stability 

O '• of the code was significantly improved as well. Thanks to these changes, PHOTOS generates 

^ [ bremsstrahlung corrections in Z and W decays with a precision of 0.1%. This precision was 

I established in El with the help of a multitude of distributions and of a specially designed test 

P,; (SDP). 

p ■ In this note we will not repeat a discussion of the design properties, but we will recall the 

O , main tests that document robustness and flexibility of the PHOTOS design. This aspect may be 

of broader use and may find extensions in future applications, also outside the simple case of 
. ^ I purely QED bremsstrahlung in decays. 

^ ■ We begin with an informal presentation of the components of the PHOTOS algorithm us- 

■ ing operator language. The consecutive approximations used in the construction of the crude 

distribution for photon generation, and the correcting weights used to construct the physically 
complete distributions are listed, but will not be defined in detail, because of limited space 
allowed for the workshop contributions. Instead, we present the variations of the algorithm. 
Comparisons between different options of the algorithm provide an important class of technical 
tests, and also helps to explore limits of the universality of the PHOTOS solution. The results of 
some of these tests will be listed later in the contribution (for the remaining ones and details we 
address the reader to refs. El SI). In the comparisons we use the SDP universal test based on 
MC-TESTER Q as in ref. Q- We must skip the repetition of its definition here as well. 

The starting point for the development of PHOTOS was the observation that, at first order, 
the bremsstrahlung corrections in the Z /J^/J^ process can be written as a convolution of the 
Bom-level distribution with the single-photon emission kernels for the emission from /j+ and 
fj . 
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The fomiulae for the emission kernels are 3 -dimensional and can be parametrized using the 
angles and the invariant mass, which are the same variables as those used in the parametrization 
of the three-body phase space (the kernels use only a subset of the complete set of phase-space 
parametrization variables). The remaining two angular variables in the kernels can be identified 
as the angles defining the orientation of the /j+ and/or /j directions (for a detailed definition, 
see e.g. El). 

The principle of the action of the single-photon algorithm working on n-body decay is to 
replace a point in the n-body phase space ^^2, with either the point in the original Q.2, or the 
point in the {n + l)-body phase space (with generated photon). The overall normalization 
of the decay rate has to change as well and, for example, in the case of Z ^ A'^A' it needs to be 
multiplied by a factor of 1 + 

Subsequent steps of the PHOTOS algorithm are described in terms of the evolution oper- 
ators. Let us stress the relations of these operators to the matrix elements and phase-space 
parametrizations. We will present the decomposition of the operators in the top-down order, 
starting with the definition of Ra, the operator describing the complete PHOTOS algorithm for 
single emission (which at least in the case of Z and leptonic x decays originates from field 
theory calculations without any approximation). Then, we will gradually decompose the oper- 
ators (they differ from decay channel to decay channel) so that we will end up with the single 
well-defined, elementary operator for the emission from a single charged particle in the final 
state. By aggregation of these elementary operators, the Ra may be reconstructed for any decay 
channel. Let us point out that the expression of theoretical calculations in the form of operators 
is particularly suitable in computer programs implementation. 

We cannot present here a separate discussion of the factorization properties, in particular 
to define/optimize the way the iteration of R's is performed in PHOTOS. Not only the first-order 
calculations are needed, but also higher-order ones, including mixed virtual-real corrections. 
For practical reasons, the Ra operator needs to be regularized with the minimum energy for 
the explicitly generated photons: the part of the real-photon phase space, under threshold, is 
integrated, and the resulting factor is summed with the virtual correction. 

• 1 

Let us define the five steps in Ra separation. In the first one, the Ra is replaced by (we 
use two-body decay as an example) Ra = Ri{Rs{lJ^) +Rs{l-i )), where Rj is a generalized 
interference operator, and Rs is a generalized operator responsible for photon generation from 
a single, charged decay-product. Here we understand the generalized interference operator 
as shifting between different kinematical configurations while respecting energy-momentum 
conservation, thus also overall normalization of the distribution under construction. 

There is a freedom of choice in the separation of Ra into Rj and Rs- The Rs operator acts 
on the points from the phase space, and the results of its action belong either to ^^2 or to 
The domain of the Rj operator has to be ^^2 + ^^35 and the results are in either Q.2 or Q3. In 
our solution we required that Rj acts as a unit operator on the t22-part of its domain and, with 
some probability, returns the points from back to the original points in Q.2, thus reverting 
the action of the Rs- 

Let us stress that in practical applications, to ease the extension of the algorithm to "any" 
decay mode, we used in PHOTOS a simplification for Rj . Obviously, the exact representation 
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of the first-order result would require Rj to be decay-channel-dependent. Instead, we used 
an approximation that ensures the proper behaviour of the photon distribution in the soft limit. 
Certain deficiencies at the hard-photon limit of the phase space appear as a consequence, and are 
the subject of studies that need to be performed individually for every decay channel of interest. 
The comparisons with matrix-element formulae, as in |l6|, or experimental data, have to be 
performed for the sake of precision; they may result in dedicated weights to be incorporated into 
PHOTOS. In principle, there is no problem to install a particular decay-channel matrix element, 
but there has not been much need for this yet. So far, the precision of the PHOTOS algorithm 
could always be raised to a satisfactory level by implementing some excluded parts of formulae, 
being the case of W decay ^ an exception. 

The density generated by the Rs operator is typically twice that of real photons all over the 
phase space; it can also overpopulate only those regions of phase space where it is necessary 
for Rj. The excess of these photons is then reduced by Monte Carlo with the action of Rj. 

• 2 

In the next step of the algorithm construction, we have separated Rs = RbRa, where Rb 
was responsible for the implementation of the spin-dependent part of the emission, and the Ra 
part was independent of the spin of the emitting final-state particle. Note that this step of the 
algorithm can be performed at the earlier stage of generation as well, that is before the full 
angular construction of the event. Rb is again, as Rj, of the generalized interference type; it 
moves the hard bremsstrahlung events in excess back to the original no-bremsstrahlung ones. 
Rb operates on the internal variables of PHOTOS rather than on the fully constructed events. 

• 3 

The definition of the Rj, Rb, Ra operators was initially based on the inspection of the first- 
order matrix elements for the two-body decays. In the general solution for Ra, the process of 
multiple-body decay of particle X is temporarily replaced by the two-body decay X chY, in 
which particle X decays to the charged particle ch, which "emits" the photon, and the "spec- 
tator system" Y. The action of the operator is repeated for each charged decay product: the 
subsequent charged particle takes the role of the photon emitter ch; all the others, including the 
photons generated in the previous steps, become a part of the spectator system Y. The indepen- 
dence of the emissions from each charged product then has to be ensured. This organization 
works well and can be understood with the help of the exact parametrization of multibody phase 
space. It is helpful for iteration in multiple-photon emission. It also helps to implement some 
genuine second-order matrix elements. This conclusion can be drawn from an inspection of the 
second-order matrix elements, as in Q. 

• 4 

In the next step, we decompose the Ra operator, splitting it in two parts: Ra = RqRx- The 
Rx operator generates the energy of the (to be generated) photon, and Ra generates its explicit 
kinematical configuration. 

The Rx operator acts on points from the phase space, and generates a single real number 
x; the Ra operator transforms this point from ^2 and the number x to a point in ^^3, or leaves 
the original point in Qi- Note that again, as Ri, the Ra operator has to be unitary and has to 
conserve energy-momentum. The Rx operator does not fulfill these criteria. 

An analogy between Rx and the kernel for structure-function evolution should be mentioned. 
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However, there are notable differences: the x variable is associated more with the ratio of the 
invariant mass of decay products of X, photon excluded, and the mass of X, than with the 
fraction of energy taken away by the photons from the outgoing charged product ch. Also, 
can be simplified by moving its parts to Ra, Rs or even Rj. Note that in Rj^ the contributions of 
radiation from all charged final states are summed. 

• 5 

The Rx operator is iterated, in the solutions for double, triple, and quartic photon emission. 
The iterated R^ can also be shifted and grouped at the beginning of the generation, ignoring 
phase-space constraints. The iterated R^ takes a form similar to a formal solution for structure- 
function evolution, but with exceptionally simple kernels. The phase- space constraints can then 
be introduced later, with the action of the Ra operators. Because of this, the iteration of Rx can 
go up even to infinite order. The algorithm is then organized in two steps. At first, a crude 
distribution for the number of photon candidates is generated; then, their energies are defined. 
At this stage we can perform a further separation: Rx = RfRo, where the Rq operator determines 
whether a photon candidate has to be generated at all, and Rf defines the fraction of its energy 
(still without energy-momentum-conservation constraint). From the iteration of Rq, we obtain a 
Poissonian distribution, but any other analytically solvable distribution would be equally good. 

The overall factor, such as 1 + | ^ in Z leptonic partial width, does not need to be lost. It 
needs to find its way to the Rq, and affects the total rate of the process. For the case of the FSR, 
discussed here, we can skip this point; however, it may be important for generalizations. 

The input data for the algorithm are taken from the event record, the kinematical configura- 
tions of all particles, and the mother-daughter relations between particles in the decay process 
(which could be a part of the decay cascade) should be available in a coherent way. 

This wraps up, a basic, presentation of the steps performed by the PHOTOS algorithm. For 
more details see dlHl- 

Tests performed on the algorithm: 

1. The comparison of PHOTOS running in the quartic-photon emission mode and the expo- 
nentiated mode for the leptonic Z and W decays may be found on our web page which 
documents the results of the tests |4J. The agreement in branching ratios and shapes of 
the distributions is better than 0.07% for all the cases that were tested. It can be con- 
cluded that changing the relative order for the iterated Rq and the rest of Ra operators 
does not lead to significant differences. This test, if understood as a technical test, is 
slightly biased by the uncontrolled higher-than-fourth-order terms which are missing in 
the quartic-emission option of PHOTOS. Also, the technical bias, due to the minimal pho- 
ton energy in generation, present in the fixed-order options of PHOTOS may contribute to 
the residual difference. 

2. The comparison of PHOTOS with different options for the relative separation between Rj 
and Rs- The tests performed for the fixed-order and exponentiated modes indicated that 
the differences in results produced by the two variants of the algorithm are below the level 
of statistical error for the runs of 10^ events. In the code these two options are marked 
respectively as VARIANT-A and VARIANT-B. 
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3. The comparisons of PHOTOS with different algorithms for the implementation of the Rj 
operator. In PHOTOS up to version 2.12, the calculations were performed using internal 
variables in the angular parametrization. This algorithm was limited to the cases of de- 
cays of a neutral particle into two charged particles. In later versions, the calculations are 
performed using the 4-momenta of particles, hence for any decay mode. The tests per- 
formed for leptonic Z decays indicated that the differences are below the statistical error 
of the runs of 10^ events. 

4. Comparison of PHOTOS with different options for the relative separation between Rq and 
Rx, consisting of an increase in the crude probability of hard emission at Rq. The tests 
performed for the exponentiated mode of PHOTOS indicated that the differences are below 
the statistical error of the runs of up to 10^ events. 

5. The remaining tests, including new tests for the effects of the interference weights in 
cascade decays, are more about the physics content of the program than on the technical 
or algorithmic aspects. They are presented in ref. lU and the results are collected on the 
web page i4|. 

Multiple options for PHOTOS running and technical compatibility of results even for 10^ 
event samples generated in a short CPU cycle time are encouraging. They indicate the potential 
for algorithm extensions. Note that PHOTOS was found to work for decays of up to 10 charged 
particles in the final state. 
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